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In this paper three things are done, (i) We investigate the analogues of Cerenkov radiation for 
the decay of a superluminal neutrino and calculate the Cerenkov angles for the emission of a photon 
through a W loop, and for a collinear electron-positron pair, assuming the tachyonic dispersion 
relation for the superluminal neutrino. The decay rate of a freely propagating neutrino is found to 
depend on the shape of the assumed dispersion relation, and is found to decrease with decreasing 
tachyonic mass of the neutrino, (ii) We discuss a few properties of the tachyonic Dirac equation 
(symmetries and plane- wave solutions), which may be relevant for the description of superluminal 
neutrinos seen by the OPERA experiment, and discuss the calculation of the tachyonic propagator, 
(iii) In the absence of a commonly accepted tachyonic field theory, and in view of an apparent 
"running" of the observed neutrino mass with the energy, we write down a model Lagrangian, 
which describes a Yukawa-type interaction of a neutrino coupling to a scalar background field via a 
scalar-minus-pseudoscalar interaction. This constitutes an extension of the standard model. If the 
interaction is strong, then it leads to a substantial renormalization-group "running" of the neutrino 
mass and could potentially explain the experimental observations. 

PACS numbers: 95.85.Ry, 11.15.-q, 03.70. +k, 05.10.Cc, 03.65.Pm 

I. INTRODUCTION 

For subluminal particles ("tardyons"), the dispersion relation for the energy E in terms of the velocity v is given 
by E = m/y/l — v 2 (with v < 1), and for superluminal particles ( "tachyons" ) , it reads as E — ra/^/v 2 — 1 with v > 1. 
Therefore, the "light barrier" at v = 1 (we set the speed of light equal to one) looks like an (infinitely) elevated 
mountain in terms of the energy of a relativistic particle. Recami pQ quotes Sudarshan with reference to an imaginary 
demographer who studies population patterns on the Indian subcontinent: "Suppose a demographer calmly asserts 
that there are no people North of the Himalayas, since none could climb over the mountain ranges! That would be 
an absurd conclusion. People of central Asia are born there and live there: they did not have to be born in India and 
cross the mountain range. So with faster-than-light particles." 

In the early morning hours of 23 February 1987 (at 2'52'36"), an unexpected neutrino bunch arrived at the LSD 
detector under the Mont Blanc roughly 4.5 hours before the rest of the neutrinos from SN1987A, and before the 
supernova became visible [5] . We are currently facing mounting evidence that neutrinos may be genuinely superluminal 
particles ("tachyons"). The MINOS experiment [3] has measured superluminal neutrino propagation velocities which 
differ from the speed of light by a relative factor of (5.1 ± 2.9) x 10 -5 at an energy of about E v « 3 GeV, supporting an 
earlier FERMILAB experiment where the trend of the data also pointed toward superluminal neutrinos [4] . This result 
has recently been confirmed by OPERA [5 with better statistics and in a wider energy interval, as detailed below. 
One of the prime candidates for a genuinely superluminal particle is the neutrino, which has never been observed at 
rest. A number of experimental groups have measured negative mass squares for the electron neutrino from tritium 
beta decay endpoints [MS] with mean values in the interval — 147 eV 2 < m 2 < for the electron neutrino mass square, 
at an energy of the order of E v « 18keV. While some recent measurements indicate values consistent with a vanishing 
neutrino mass [TUHT2"] at even lower energies, the mean value of the experimental data for m 2 (electron neutrino) still 
is negative and of the order of a few negative eV 2 (for an excellent overview, see Ref. [IS])- The idea that neutrinos 
might be of tachyonic character is not new [i"4Tfl7] . Tachyonic neutrinos fulfill the dispersion relation E 2 ~p 2 = —m 2 
with an (initially) constant parameter m u . The quantity — m 2 can be interpreted as the negative mass square of the 
neutrino. The current situation indicates the need for a convenient descriptions of tachyonic fermions. 

Ever since the early days of relativity, the notion of superluminal propagation has intrigued physicists |18) , and the 
name "tachyon" was eventually coined in Ref. |19j . The main problem in the description of a quantum field theory 
with superluminal propagation is not the superluminal velocity itself [20] , but the construction of field operators and 
the time ordering, which is in disarray because the time ordering of two space-time points which are separated by a 
space-like interval is not invariant under (subluminal) Lorentz boosts. Generally, it has been assumed that any particle 
in relativistic quantum theory should be described by a unitary irreducible representation of the Poincare algebra or 
its supersymmetric generalization. It may be necessary to relax this restriction somewhat in order to accommodate 
for a field theory of supersymmetric tachyons [2"TH2"3"] . Three recent review articles [TJ [2U [25] provide rather detailed 
background information on the development of the theory of superluminal particles. 

The recent OPERA experiment [3] uses a baseline of L = (731278.0 ± 0.2) m. Two clocks used in the measurement 
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are accurately synchronized by a technique used to compare atomic clocks [35] [37] . It is of particular importance that 
the synchronization of the two systems was calibrated by the Federal Swiss Metrology Institute METAS (Bundesamt 
fiir Metrologie) in 2008 and verified in 2011 by the Federal German Metrology Institute PTB (Physikalisch-Technische 
Bundesanstalt). As reported in Ref. jS], the difference between the time base of the CERN and OPERA receivers was 
measured to be (2.3 ± 0.9) ns and is taken into account in the evaluation of the measurement. The four data bins are 

E v = 13.8 GeV , St = (54.7 ± 18.4 ± 7.1) ns , A = ^— C = (2.24 ± 0.75 ± 0.29) x 10~ 5 , (la) 

c 

E v = 28.2 GeV , St = (61.1 ± 13.2 ± 7.1) ns , A = ^— C = (2.50 ± 0.54 ± 0.29) x 10~ 5 , (lb) 

c 

E v = 40.7 GeV , St = (68.1 ± 19.1 ± 7.1) ns, A = = (2.53 ± 0.78 ± 0.29) x 10" 5 , (lc) 

c 

and the overall average is 

E u = 17 GeV, 5t= (57.8 ±7.2 ±7.1) ns, A = = (2.37 ± 0.32 ± 0.29) x 10~ 5 . (Id) 

c 

While the OPERA data rather point to a slight increase in the ratio A = (v ~ c)/c with the neutrino energy, than to 
a trend in the opposite direction, the data are generally consistent with a constant ratio A = (v — c)/c in the entire 
energy interval 13.8 GeV < E v < 40.7 GeV. 

Tachyonic neutrinos fulfill the space-like dispersion relation E 2 — p 2 = —m 2 and travel faster than light. Super- 
luminality is conserved under Lorentz boosts (see Ref. [20] and Fig. [2] below) . It has been argued that neutrinos 
traveling at velocities consistent with the recent OPERA data should decay by neutral massive analogues of Cerenkov 
radiation [28] . The noncovariant dispersion relation E v — \p\ v v has been used in recent work on the subject |28j (here, 
v „ denotes the neutrino velocity) . Freely propagating subluminal relativistic particles as well as tachyons [T] [2U [2"S"] 
fulfill the "opposite" relation |p| = E v v v . Both relations E v — \p\v v and \p\ = E v v v lead to a large virtuality 



\E 2 — p 2 \ on the order of (117 MeV) 2 when applied to the recently measured OPERA data [see Eqs. (12) and ( fl3| ) 
below] . These observations are inconsistent with beta decay end point measurements [STU2| which have led to values 
of a few eV 2 , for neutrinos in the keV energy range. This confusing situation raises a number of questions. Starting 
from the tachyonic Dirac equation, we conclude that additional interactions, hitherto not accounted for, arc required 
in order to explain the OPERA data which exhibit a larger-than-expected virtuality at higher energies, or, expressed 
differently, an energy-dependent mass. 

At the current, early stage in the development of theoretical models describing superluminal particles, a certain 
degree of speculation cannot be avoided. For completeness, we should note that we neither consider models based 
on deformed special relativity [2"9"H3"2"] nor kinematic constraints resulting from such models [28;, 33"H3"5] in any greater 
detail. Lorentz symmetry is conserved in our approach. 

We start with a digression on the kinematic constraints to the observation of neutrinos along the OPERA baseline 
in Sec. In] The tachyonic Dirac equation and its solutions are being reviewed in Sec. |III[ Chiral Yukawa interactions, 



which induce a neutrino mass running via the renormalization group (RG), are studied in Sec. IV Conclusions are 
reserved for Sec. |vj We always carefully distinguish between \p\ and the four-vector p, and we use natural units with 
H = c = eo = 1. 



II. KINEMATIC CONSTRAINTS 



The recent OPERA experiment has analyzed the propagation of muon neutrinos. If neutrinos propagate faster 
than the speed of light, then a number of decay processes are kinematically allowed which are otherwise forbidden. 
These include the following decays (see Fig. [I]), 

"jt -> + 1 , (2a) 
-> V,,, + e + + er , (2b) 
->• Vf,, + v e + v e . (2c) 

In Ref. |28j . these decay processes are analyzed under the assumption of the Lorentz- violating dispersion relation 

dE 

-— - = const. , E v = \f„\ v v , v v rj 1 + A , (3) 



(a) (b) (c) 

FIG. 1: Feynman diagrams for the decay processes of a tachyonic superluminal neutrino, as given in Eq. (pj). The tachyonic 
neutrino may emit of photon via a W loop [Fig. (a)], or an electron-positron pair, [Fig. (b)], or a neutrino-antineutrino pair 
[Fig. (c)]. The processes scale with the quantum electrodynamic (QED) coupling constant a and the weak coupling constant 
Gf as follows, (a) is proportional to aG%, (b) is proportional to oGf, and (c) is proportional to a 2 Gp. 

where A = 2.37 x 10~ 5 corresponds to the value given in Ref. [5]- Processes (a) and (c) are parametrically suppressed 
with respect to process (b), and therefore process (b) is deemed to constitute the dominant decay channel. 

One may observe that the dispersion relation E v = \p v \ v v is at variance with both the subluminal (also called 
tardyonic, see Ref. |19j ) dispersion relation for freely propagating massive neutrinos, 

m„ . . mv v 

E v = r — ==, \p v \ = —=== E v v v , v u <l, (4) 

V 1 — v v V*-~ V v 

as well as with the dispersion relation for superluminal (tachyonic) particles [U Ell [25, 36 38], which reads 

Ev = To \ ' %\~ To = \ = E V V V , V V >1. (5) 

v v i - 1 v v i - 1 

In both cases Q and ([5]), one obtains \p v \ = E v v v , not the opposite relation E v = \p v \v v used in Ref. [55]. Under 
Lorentz transformations, superluminality of tachyonic particles is conserved (see Fig.[2|. In two recent papers [3^114*0] . 
it has been observed that the conclusions of [28 would change if the dispersion relation were different. Here, we are 
concerned with a more general question: Namely, to investigate how the kinematic constraints change when we 
assume a tachyonic dispersion relation for the neutrino, and whether the process ([2| is still kinematically allowed 
when El — p*£ < 0. 



For the process (2a), an easy calculation based on the energy and momentum conservation conditions reveals that 
E V =E' V + E 1 , p v =tf v + k 1 , E v = _ m 2 j e' v = a/pL 2 - ml , E 7 = \k 7 \ . (6) 
Squaring the energy conservation condition, one obtains 

El = ($ v + fe 7 ) - ml = p" u 2 + fc 7 2 - ml + 20 u ■ % , (7a) 
El = (E' u + E 7 ) 2 = ^ + fc 7 2 -ml + 2 |fc 7 | VpI 2 - ml , (7b) 

jf v 'ky = |fc 7 | ^JpJ- - ml . (7c) 
We conclude that under the assumption of the Lorentz-covariant, tachyonic dispersion relation ([5]), weak-interaction 



Cerenkov radiation is allowed. In view of Eq. (7c), the photon is radiated off at a Cerenkov angle 



C0S e 7 = i^= ^i- mi = §- = l F< i, (8 ) 
7 |fc 7 ||#J lft I IpII K 



under the assumption of a tachyonic neutrino with dispersion (7a). One may add that the kinematic consideration 
is somewhat analogous to that for the emission of ordinary Cerenkov radiation. The important observation is that 
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FIG. 2: (Color online.) Illustration of the Einstein velocity addition theorem w — (u + v)/(l + u v), in the superluminal domain 
with u £ (—1, 1) and v £ (1, 3). For superluminal v, the range w 6 [—1, 1] of values is excluded, as shown by the rectangular 
box. 



under the tachyonic dispersion relation (5j, the emission of a photon by the neutrino is always allowed, i.e., there is no 
threshold energy for the neutrino and there is no threshold for the tachyonic mass — m 2 . Once the particle becomes 
tachyonic, weak Cerenkov radiation is kinematically allowed, but the Cerenkov cone narrows as — m 2 —> 0. For a 
particle fulfilling the noncovariant dispersion relation E v = \f? v \v Vl with v v > 1, the modified Cerenkov angle cosO'^ 
is easily computed as 

assuming a neutrino with the dispersion E' v = p' v v' v and v' v > 1. This is very well approximated by cos#^ f=a l/v[, for 
v' v « 1. 

As a second step, let us consider a process in which a tachyonic neutrino fulfilling Eq. ^ emits a massive neutral 
vector meson of mass mo- This is not depicted in Fig. ([I]) but still instructive. The kinematic conditions change, 



E V = E' V +E , $ v =$ v + ko, E v = - ml , E' v = ^tfj 2 - ml , B = V4 2 + m o- ( 10 ) 
The Cerenkov angle then becomes 

ml + 2 Jfeg + ml ^Jffi - ml J fcg + -y/^ 2 - ml 

cos^n = ~ =; > (11) 

2|fco|KI IfcolKI 

where the last expression is valid in the high-energy limit, i.e, for |fc | 3> m , and 3> m^. If the vector meson 
carries away the bulk of the energy, i.e. Eq — xE v and E' v = (1 — x) E u , with x < 1, then for highly energetic incoming 
superluminal neutrinos, one can always find a narrow cone near #o in which vector meson emission is possible. 
Again, for highly energetic tachyonic superluminal neutrinos, we conclude that there is no kinematic constraint on the 
size of the tachyonic mass term —m\ which would restrict massive vector meson emission. Once the particle becomes 
tachyonic and the energy of the tachyonic particle is large enough, massive vector emission becomes kinematically 



allowed in a narrow angular region. By contrast, if we replace in Eq. (11) — mf, — ¥ +m^, we would have cos#o > lj 
forbidding vector meson emission. Also, the Cerenkov angle 6q vanishes in the limit m„ — > 0. Using more extensive 
calculations, we have checked that the same statement applies to the light fermion pair emission given in Eq. (2b) 
and depicted in Fig. [ijb). Ccrcnkov-type pair emission becomes kinematically possible for highly energetic neutrinos, 
in a narrow angular region. 

In the application of the tachyonic dispersion relation ^ to the OPERA data, we face a dilemma which also plagues 
the application of the Lorentz-noncovariant dispersion relation Namely, we have for the OPERA data according 



to Eq. (Id) 



il = E 2 V - vl = El [1 - (1 + A) 2 ] = - (117MeV) 2 [dispersion relation §] (12) 
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which is at least six orders of magnitude larger than the neutrino masses at low energy [BT-112]. Likewise, assuming 
the dispersion relation ^ implies that 

El - vl = vl [(1 + A) 2 -l]&El [(1 + A) 2 - 1] = (117 MeV) 2 [dispersion relation ^} . (13) 

If we define the expression |_E 2 — p 2 | as the "virtuality" of the neutrino which measures the deviation of the neutrino 
propagation velocity from the speed of light, then we can say that at the high OPERA energies, the neutrino velocity 
was not expected to deviate so much from the speed of light, neither in the superluminal nor in the subluminal 
direction. For example, if OPERA had hypothetically found a result of 



A = —2.37 x 10 5 [opposite sign as compared to Eq. (Id)] , (14) 



then this would have been equally surprising. In the latter case, one would probably have concluded immediately 
that the neutrino must be subject to a hitherto unknown interaction at high energy, modifying its effective (running) 
mass. We here advocate the viewpoint that the same conclusion should be drawn from the OPERA data: namely, the 
neutrino is genuinely tachyonic and subject to an unknown interaction at high energy which modifies its mass and its 
decay channels. Otherwise, it seems that the high-energy OPERA data [5] (in the GeV range) cannot be reconciled 
with the low-energy experimental results (in the keV range) of Refs. [STU2j. Of course, this statement holds provided 
the OPERA data are not subject to a hitherto undiscovered systematic error. 

The data bins given in Eq. are consistent with an energy-independent propagation velocity. While the quantity 
A need not be energy independent over large energy intervals, it appears to be so in in the energy interval 13.8 GeV < 
E v < 40.7 GeV. Therefore, in this energy interval observed by OPERA [5], the dispersion relation is assumed to be 
close to a linear relationship 

m v =m(E v )KnE v , 13.8 GeV < E v < 40.7 GeV , n = J (I + A) 2 - 1 = 6.88 x 10~ 3 « • (15) 

145 

The unknown interaction leading to the energy-dependent mass must now be investigated. When calculating de- 
cay rates, the existence of the additional interaction implies that one should use eigenstates of the neutrino in the 
additional, hitherto unknown interaction potential (i.e., taking into account the running mass) rather than freely 
propagating tachyonic states, with an effective energy-dependent tachyonic mass m v = m v (E u ) oc E u . 



III. TACHYONIC DIRAC EQUATION 

Given the obvious inconsistency of the OPERA data [5] with low-energy neutrino data [BHU] > as manifest in the 



energy-dependent effective mass (15), one may ask why an equation that describes a genuinely tachyonic neutrino 
with an energy-independent, fixed tachyonic mass m„ should be considered at all in the following. The reason is that 
if the neutrino is genuinely tachyonic, then one has to start from an equation which describes a genuinely tachyonic 
particle, with the possibility to describe additional perturbative interactions that modify the high-energy behaviour. 
Expressed differently, we would expect the tachyonic Dirac equation given below to describe the low-energy behaviour 
of neutrinos [6T-TT2"] , while the large deviation from the light cone seen at high energies [31 [5] should be ascribed to 
additional interactions. We briefly recall here that the Lorentz-covariant tachyonic Dirac equation reads 

(i7»<wm,)«*>=o, 7"=(Y-iL). '-(-^S)' 7S =UVr)' < 16 > 

Here, x = (t, x), and we use the Dirac matrices in the Dirac representation [T7]- The partial derivatives are = d/dx 11 , 
while 7 5 = i7° 7 1 7 2 7 s is the fifth current matrix. The tachyonic Dirac equation has been briefly discussed in Ref. [T4T - 
[16] . It has recently been verified that this equation is CV, as well as T invariant [17]- These symmetry properties 
apply to neutrinos. The positive-energy plane-wave solutions [17] of the tachyonic Dirac equation have the properties 

*(a;) = -$=U ± (k v )e- ik "- x , k u = (E v ,k u ), E v = Jk*-ml, \k v \>m v . (17) 



The negative-energy solutions |17j are given by 



= -$=V±(k v )j k "- x , K = {E v ,k v ), E v = Jkl-ml, \K\>m v , (18) 



V 



() 



where V is the normalization volume. These states are normalized, with Ut(k u ) U+(k„) = U^(k v )U-(k v ) 
V£(k v ) V+(k v ) = V*(k v ) V-(k v ) = 1. The spinors entering these expressions read as 



/ m„ - E u + I k v I 



i/2y/{E v -\kv\y + ml 
m u + E„ — \k u \ 



\ y/2yJ{E v -\ku\) 2 - 



( m v + E v - \kJ 



MM = 



J 



V2^(E u -\k„\) 2 + mi 
—m v + E v - \k„ 



a-{k v ) 



\ 



a—{k v ) 



(19) 



V V2^{E v -\k v \Y+ml ) 



where the helicity spinors a±{k v ) are given below in Eq. (211. If we are interested in the massless limit, then we 
should first take into account the fact that massless particles propagate at velocities very close to the light cone. For 
v = 1 + A, we have E — \k u \ ~ — m A/2 <C m. Therefore, letting A — ► 0, the dominant term for the massless limit 

a+{k v ) 

negative-energy eigenstates are given by 



actually is the mass m 3> E — \k v \. This implies, e.g. 



that U+(k v ) -> ^ ( 



for the massless case. The 



/ -m v - E v + \k v \ 



MM 



V2j(E u -\k„\) 2 + ml 



a + (k v ) 



\ 



\ s/2J(E v -\k v \Y 



a+(M 



( -m v + E v - \ie v \ 



MM 



J 



V2^/(E u -\k v \) 2 + m 2 
\L\ 



a-(k v ) 



\ 



= a-(k v ) 

\V2j(E v -\t\) 2 + m 2 J 



(20) 



The helicity spinors entering these expressions are given in terms of the polar and azimuthal angles 9 and ip of the 
three- vector fc„, 



a+(k v ) 



(I) 



sin I 



a_(M 



sin (|) e- iv 
cos ( |) 



(21) 



and fulfill 



a ■ k v 



a+(k v ) = a + (M , 



a ■ k v 

m 



a_(M = —a + (k„) ■ 



(22) 



For plane waves, E v = yk 2 — m 2 , and p v = k v fulfill the tachyonic dispersion relation which we recall for 
convenience, 



E v v v , 



v v >l, 



(23) 



so that yjk 2 - m 2 , never becomes imaginary. For k 2 < m 2 , we have resonance and antiresonance energies. We start 
with the resonances, whose energies have a negative imaginary part, 



MM = 



E„ 



( m„ + ±r„ + \k v \ . r , 

j MM 

V2 Jk 2 + m 2 + \Tl 



\ 



a+{k v ) 



MM 



( m v - hT v - \k v \ r \ 

, = MM 

V2Jk 2 + m 2 + lT 2 



-m v - kT v - | 



V V2Jk 2 +m 2 + \Y 2 



a_(fc„) 



M = -i\Jm 2 v -kl, 



k 2 < ml . 



(24a) 



(24b) 
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The antiresonance energies have a positive imaginary part, 



S+(k v ) 



( -m v - ir v + \k v \ \ 
— — -, a+\K v ) 

/ 



2 1 " 



V V2 J*a + m 2 + | 



S-{k v ) 



V2^/k 2 + m 2 



1 p2 

!/ 1 4 - 1 - w 



0_(fc«) 



o_(fc J/ ) 



= h r„ — i 



1.2 



fe,? < ml 



(25a) 



(25b) 



These states are also normalized, with R+(k v ) R+(k v ) = Rt(k v ) R-{k v ) = S+{k y ) S + {k u ) = S~t{k u ) S-(k v ) = 1. 
The term "resonances" is used in the physics literature in two contexts: (i) in order to designate the complex energy 
eigenvalue of a Hamiltonian, and (ii) in order to designate the peak in a cross section or a quantum state which 
can decay into a final state with a different particle content. In the current case, the interpretation (i) is relevant. 
The resonances have complex resonance energies; the waves are evanescent (exponentially damped) just like the 
diffracted wave under total reflection, or a wave in a waveguide below the minimum frequency for the TEi o mode 
necessary for propagation, and the resonance energies are complex just as in the case of a resonance energy of the 
Stark effect |4I) . Resonances are damped for propagation forward in time, antiresonances for propagation backward 
in time, in accordance with the Feynman prescription. The wavelength of the resonance states is too long to be 
supported in a genuinely superluminal wave packet of tachyonic mass m 2 . 
The noncovariant, Hamiltonian form of Eq. (16 1 reads as 



Hsijj(x) = (a • p + P 7 5 m„) ip(x) = E v tp(x) , 
where /3 = 7 , and a — 7 7. The Hamiltonian H§ has the pseudo-Hermitian 42- 50 property 

H = V H+{x) V- 1 = PH+{~x) P- 1 , 



(26) 



(27) 



where P is the full parity transformation and P is the parity matrix P = 7 . The eigenvalues of a pseudo-Hermitian 
operator come in complex-conjugate pairs and are real if the tachyonic dispersion relations ^ are fulfilled. This 
can be seen as follows. Because the spectrum of a Hermitian adjoint operator consists of the complex conjugate 
eigenvalues, we have an eigenvector <p(x) with eigenvalue E* provided there exists an eigenvector ip(x) with eigenvalue 
E, 



H 5 {x) ip(&) = E4>(x) , H+(x) (j)(x) = E* (j)(x) . 
Then, the transformation x — > —x and the introduction of the P = 7 parity matrix leads to 

H+{-x) <j){-x) = E* <p(-x) , PH+i-^p- 1 (P(f>(-x)) = E* P<j)(-x) 
By assumption, PH^(— a ? )_P~ 1 = H${x) and thus 



H 5 (x) P(j){-x) = E* P(j>(-x) , H 5 (x) ip{x) = E*ijj(x) 



tp(x) = P(f>(-x) 



(28) 
(29) 

(30) 



This implies that ip{x) = Pip(-x) is an eigenvector with eigenvalue E* . The eigenvalues of i?5 thus come in complex- 
conjugate pairs, and furthermore, they are real for plane waves fulfilling the dispersion relation ([5]). 
The covariant Green function corresponding to the Hamiltonian H§ thus reads as 



Sr(p) = 7° 



E-H K 



p + 7 5 m v 



(31) 



The tachyonic poles at E„ —p 2 = —m 2 , have to be encircled in a way consistent with the boundary conditions imposed 
on the Green function. Eigenvalues with E 2 = p 2 — m 2 < represent evanescent waves. If one encircles the poles of 



the Green function according to the Feynman prescription, 



S T (p) 



1 



7 5 (m v + i e) p 2 + m 2 + 



(32) 



then the energy-momentum dispersion relation is infinitesimally displaced to read E v = ±-v/p^ — mf 
consistent with the evanescent wave picture because positive-energy solutions have the form E v = e 



— i e. This is 
i\/V 2 ~ m l\ 



and are thus exponentially damped for the propagation into the future, whereas negative-energy solutions have the 
form E v = — e + i \J\p 2 — m 2 \ and are thus exponentially damped for the propagation into the past. In general, the 
Feynman prescription assigns an infinitesimal negative imaginary part to energies whose real part is positive, and vice 
versa. 

Thus, while the time propagation of strictly tachyonic wave packets (superpositions of the tachyonic plane-wave 
solutions) is fully unitary (they have real eigenvalues), a slight violation of unitarity cannot be avoided if one allows 
eigenstates with E 2 = p 2 — m 2 < 0. The complex resonance energies (the real part is only infinitesimal) 



E„ 



iv4 



P- 



nl\, E v = -e + i^\p 2 -ml\, p 2 < m 2 v , (33) 

describe the suppression of subluminal components of a superluminal wave packet under time evolution. One has to 
allow these solutions in the propagator (32 1 if one would like to carry out the Fourier transformation consistently, i.e., 



over the entire range p 11 £ R , or describe the time evolution of a general wave packet under the Green function (32 1. 
It seems that a slight violation of unitarity, relevant to the small sector p 2 < m 2 , where m„ initially is on the order 
of a few eV, is a price for the introduction of tachyonic particles [33]. Note that full unitarity cannot be preserved 
anyway in a tachyonic theory if one goes beyond tree-level amplitudes, as shown in Ref. 15 1| . The time propagation 
of wave packets in potentials with manifestly complex resonance energies has been described in Refs. [521 [53] . The 
evanescence of the subluminal neutrino wave function components, which are excluded from the real neutrino plane- 
wave eigenstates but included in the propagator, is somewhat analogous to the photon propagator, where one includes 
the so-called scalar and longitudinal photons in the photon propagator but leaves them out from the real, physical 
states of the photon field, which are composed of transverse photons. 

In Ref. |17j , the tachyonic propagator ( 32 1 is derived not by inversion of the Hamiltonian, but by a quantization of 
the tachyonic field operators. We briefly sketch the essential elements of the derivation. The field operator is written 
as 



(34) 



where E v — y k 2 — m 2 — i e is the tachyonic energy and the four- vector k v equals fc„ = {E w , k v )- Here, the b operators 
annihilate particle, whereas the d + create antiparticles. We here explicitly accept a Lorentz-covariant vacuum state, 
which transforms according to Ref. [19). The Lorentz-transformed vacuum is filled with all particle and antiparticle 
states whose energy changes sign under a Lorentz transformation (Lorentz boost), as outlined in Eqs. (5.6) and (5.7) 
of Ref. [19] . The spinors U and V are given by 



U+{k v ) 



hi, 



E v + \kJ 



2 ^fm^ y'\k v \-m l 
m,, + E,, - IfcJ 



\2. 



a + {k v 
a + {k u 



\ 



for positive energy, and by 



V+(M 



k,. 



/ 



rn,, 



E„ 



2 ^vn v 



V 2 



= a + (k v ) 



- m h 
\L\ 



V 2 



'm,. 



a + (k v ) 



k<. 



v_(M 



2 y \k„\ - m l 
m,, + E,, -\k,,\ 



V 2 



a-{k y ) 
a-{k y ) 

a-(k v ) 
a-{k v ) 



\ 



I 



for negative energy (in both cases, we assume that \k v \ > rn). The normalization conditions are given by 



U a {k v ) U a {k v ) = Ut{k v )^U a {k v ) 



Quantizing the theory according to Fermi-Dirac statistics, 

{b„(ku)MK)} = {a+(k v ),a+(k' v )} = {driK^dpQi)} = {dt(K),d+(K)} = 0, 



{b a (k,),b;(kl)} = (- ( t)(2tt) 3 -5\K-k' u )8 ap , {d a {k v ),d+(k' u )} = (-a)(2ir) 3 -6 3 (kv-K)^ P , 



E 



(35a) 



(35b) 



(35c) 

(36a) 
(36b) 
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one can easily show that 



^(-<r) U a {ku)®U a {k v ) 



7 



7 5 m 



2m 



2(-<r) v CT (M®v CT (£)7 



2m 



(37) 



The field anticommutator is 
d 3 fc„ m 



(2tt)3 t 

where £ denotes the spinor index. It follows that 
{^(x)> e (y)} 7 5 = (i0-j 5 m u )^,iA{x-y) 



(38) 



A(x - y) 



d 3 fc„ 1 

(2tt) 3 2e; 



e -ifc„-(x-3/) _ e ife„-(x-y) 



where A(x — y) is introduced as in Chap. 3 of Ref. [54]. Furthermore, Eq. (3.170) of [54 finds the generalization 



(39) 



= - ( 7 °) ss , a A(x-y) _ = ( 7 °) £f , 8\x-y) 



II' 



(40) 



In full analogy with Eq. (3.174) of Ref. [SJ and in agreement with Ref. [55], the tachyonic (T) propagator is then 
found as 



T^(x)^,(yh 5 0)=iS T (x-y) 



ill' . 



St(x - y) 



d 4 k 



(27T) 



1.2 



(41) 



which confirms Eq. (32). Indeed, the propagator obtained from the quantized theory is equal to the propagator 



obtained from the inversion of the Hamiltonian in the Lorentz-covariant formulation, as it should. The couplings of the 
neutrino involve the chirality projector (1 — 7 5 )/2, and in view of 7 5 (l± 7 5 )/2 = ±(l± 7 5 )/2, the introduction of the 7 5 
matrix in Eq. (32 1 is reabsorbed into the interaction Lagrangian. The non-unitarity is small because toJ is very small. 



For a tachyonic particle, the evanescence of non-tachyonic wave packet components is natural because its tachyonic 
components remain tachyonic upon Lorentz transformation (see Fig. |2]). The tachyonic Dirac equation provides for a 
convenient framework to describe freely propagating, superluminal, electromagnetically neutral, particles. 



IV. NEUTRINO MASS RUNNING 



In Sees. [IT] and III we have seen that a running neutrino mass (with the energy) is able to conceivably suppress 
Cerenkov-type decay processes, and the quantization of the tachyonic Dirac equation has been discussed as a convenient 
description for tachyonic spin-1/2 particles; it naturally implies the suppression of the right-handed neutrino. If current 
experimental data [5] are confirmed, then we now have to explain why the effective mass of the neutrino, which needs 
to be inserted into the tachyonic Dirac equation, changes from a few eV in the keV neutrino energy range, to a mass 
on the order of MeV in the GeV energy range. We note that neutrino mass running is usually assumed to initiate 
on the energy scales of Grand Unification (see-saw mechanism). However, the experimental data 03 G>HS] all point 
to a neutrino mass running which sets in at much lower energy scales. We assume that the mass term is genuinely 
tachyonic. 

The scenario that we would like to propose is as follows: We conjecture that the neutrino mass running is due 
to an interaction with a hitherto unknown field that modifies its effective mass with the energy. At low energy, the 
interaction with the unknown field is weak, so that the apparent neutrino mass is in the eV range, whereas at higher 
energies, the interaction becomes stronger and leads to the observed [S] large tachyonic masses. We thus assume that 
the (bulk of the) neutrino mass is created dynamically [51] . Possibly, there is some threshold region where the effective 
mass of the neutrino intersects with the mass of the field it interacts with, and this might help explain consistency 
with astrophysical data [2 j. In the following, we would like to present a semi-quantitative analysis which supports 
these conjectures. 

We investigate a scalar-minus-pseudoscalar (S — P) interaction Lagrangian of the form 



(42) 
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Here, 4>x is a scalar field operator, G is a dimensionless coupling, and the fermionic field operators for the neutrino 
are denoted as ip. The operator is of dimension 4 and therefore renormalizable; it describes a Yukawa interaction 
with a chirality projector. The complete Lagrangian of the tachyonic neutrino field, the scalar field plus the S — P 
interaction reads 

C{x) = X - [fey (fyfe) - (fyfe) y$(s)] -fe 7 6 m„fe 



■ 4> x (x) (□ + M|) 4> x {x) + G4> x (x) fa) (1 - 7 5 ) fa) . 



(43) 



At low energy, from dimensional analysis alone, the induced one-loop neutrino mass running via the renormalization 
group (RG) can be written down as 



dm v dm„ 3 2 

dm(/i) d/i 



[G^)] 2 = G| ln(/i) = ^ ln( M ) 



(44) 



where we assume a logarithmic running of the coupling constant with the scale /x. Integrating the RG evolution 
equation, 



dm,, 



G 2 



x 



dfx 
I 1 



m„(17GcV) 



dm,. 



G 2 



x 



17GcV 



m„(18kcV) 



18koV A 4 



(45) 



with m„(18keV) lOOeV (see Ref. [g]) and (17 GeV) « 117MeV (see Ref. E)), we find that an X particle of mass 
in the range Mx ~ 1.4 kcV could potentially induce a neutrino mass running from about 100 eV at 18 keV energies [6] 
to 117 MeV at energies of 17 GeV [5]. Here, we assume that G ~1 and a universal running of the electron neutrino 
mass [H| and the muon neutrino mass [3] with the energy. The difference in the observed OPERA neutrino mass [3] 
of 117 MeV with low-energy neutrino data Refs. [BTH2]. where masses in the eV range were observed, suggests that 
significant neutrino mass running has to set in at energies much below 17 GeV, so that we can safely assume that 
Mx <C 17 GeV. This finding and the interaction (42) is not described by any known particle in the standard model, 



and thus, our model constitutes a pertinent extension. However, one may object that this treatment amounts to an 
application of a one-loop running of the mass in a domain which in view of G ~ 1 clearly is nonperturbative. 

This high-energy limit could be analyzed as follows. We first recall that in the high-energy domain, where the 
effective neutrino mass is in the MeV range (see Refs. [HE]) we assume that the neutrino mass is (almost) exclusively 
generated by the strong (nonperturbative) self-interaction with the X field. It is interesting to observe that polynomial 
behaviour of RG functions in the strong- coupling domain has recently been obtained by a sophisticated analysis of 
higher-order perturbative terms, for the fi functions of <f> theories and of quantum electrodynamics [571 158j . If the 
mass of the X particle is negligible as compared to the mass of the neutrino in the high-energy domain, then the 
mass scaling must be independent of Mx, and again, from dimensional analysis alone, we may conjecture that in the 
high-energy, strong-coupling limit, 



dm,. 



d/i. 



oc G m v . 



Am,, 



KG 



d/i 



K G 



(46) 



where K is & constant of order unity. In view of Eq. (15), if we assume that G ~ 1/s/K, then 

1 



m v = m v {E v ) = rj {E v ) 



KG* 



VE,, G ■ 



K 



n 



i 

145 



(47) 



where the value of ry is chosen such as to be consistent with Eq. ( 15 ) 



We are now in the position to add some more, somewhat speculative, remarks on the experimental findings of 
Refs. [3H5]. Based on the numerical entries in Table II and Fig. 3 of Ref. [4 , one may investigate the observed 
neutrino velocities as a function of the propagation energy. The authors of the somewhat inconclusive 1979 paper 
(Ref. [J]) suggest to ascribe a path length correction of A pat h = —0.5~^ 01 x 10" 4 to their data, because the muons 
that were "racing" against the neutrinos in the experiment were assumed to be artificially delayed due to multiple 
scattering events, which extend the muon path length in comparison to the muon neutrino path length. The path 
length correction was assumed to be constant over the energy range analyzed in Ref. |3] , uniformly affecting neutrinos 
in the energy range of 32 GeV < E v < 195 GeV in an experiment over a relatively short baseline of about 900 m 
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FIG. 3: (Color online.) Measured neutrino velocities in the range E v = 3 GeV (Ref. [3]) up to JS„ = 195 GeV (Ref. [4). The 
OPERA data are given in Eq. and correspond to the data bins at E = 13.8 GeV, E = 28.2 GeV, and E = 40.7 GeV (circles). 
The data point at E v — 3 GeV is from Ref. [3] (square). All remaining data points (triangles) are from Ref. [4]. Panel (a) 
corresponds to the data plotted in Fig. 3 of Ref. [3], while panel (b) applies a path length correction of A pat h = —0.5 x 10 -4 to 
the data (triangles) of Ref. [3], as discussed near the end of Ref. [3|. Here, A is the relative deviation from the speed of light in 
vacuum, which we multiply by a scaling factor 10 4 on y axis. The solid line at A = 2.4 x 10 -4 corresponds to the result 1 48 1 
based on our model. 



(which is smaller than the OPERA baseline by a factor of roughly 10 3 ). We find that the discussion on the derivation 
of the path length correction in Ref. g] is rather short and therefore present data with and without this correction 
in Fig. [3j the same approach was recently taken in Figs. 1 and 2 of Ref. [31]. The model (47) leads to a constant 
deviation of the neutrino velocity of 



A 



1 = 2.4 x 10" 



(48) 



independent of the neutrino energy. This result is compared to the available experimental data [3H5] in Fig. [3j While 
our model is somewhat speculative at the current stage, it is intriguing to observe that the solution of the simple- 
minded RG equation (47) is in good agreement with the observed neutrino velocities over a wide energy interval 
(see Fig. [3]). We also recall that the concomitant significant neutrino mass running will suppress decays because the 
tachyonic mass in the exit channel is much lower than in the incoming channel (see Sec. pl|. 



CONCLUSIONS 



Tachyons have a potential of fundamentally altering our view of physical law, but they can be incorporated into the 
framework of Lorentz transformations, despite obvious problems with the causality principle. In Ref. |36) . the authors 
argue that a "sensible" theory is obtained if one insists that the only physical quantities are transition amplitudes, 
and a negative-energy in (out) state is understood to be a positive-energy out (in) state. This statement is in need 
of further explanation. Suppose that observer A sees event E before E' , and observer A' sees event E' before E, 
because the two events are separated by a space-like interval, and the Lorentz transform for the frames A and A' 
reverses the time-ordering of events E and E' . According to Ref. |20j . the reversed time ordering occurs if and only if 
the energy between the two frames also changes sign. So, provided one reinterprets the negative-energy eigenstates of 
tachyonic Dirac Hamiltonian propagating backward in time (the antiresonances included) as positive-energy solutions 
propagating forward in time, the creation and absorption of a particle can be consistently reinterpreted if only the 
transition amplitude is unaffected by the reinterpretation. This point has also been stressed in Refs. [2TJ1 155] . 

One problem, though, in the consistency of observations of tachyons lies in conceivable decay processes [28] . In this 
paper, we investigate threshold conditions for the emission (see Sec. [TTJ) of real particles by analogues of Cerenkov 
radiation emitted by superluminal, tachyonic neutrinos that fulfill the dispersion relation (15]). We find that such 



emissions, as shown in Fig. [T 



angles 9 [see Eqs. (|8| and (11 



are possible at high energies for small Cerenkov angles in a narrow cone of emission 
Furthermore, at sufficiently large energy, a nonvanishing emission probability exists 
for even very small tachyonic mass squares — rr? v . However, the calculation of the corresponding decay rates crucially 
depends on the dispersion relation used in the calculation. The tachyonic relation ([5| is Lorentz-invariant, and the 
effective mass m„ crucially influences the decay rate. We then investigate, based on the tachyonic Dirac equation (see 
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Sec. Ill), how the effective neutrino mass m u could possibly change from a few eV at low energies in the keV range 
to energies of a hundred MeV in the GeV range [see Eqs. (12), (13) and (15)]. 

We here come to the conclusion that a viable explanation for the large virtuality E 2 — p 2 of the OPERA neutrinos 
could be due to an additional interaction that modifies the neutrino propagation at high energies. At an energy in the 
GeV range, as measured by OPERA, the propagation velocity of a particle with a rest mass on the order of a few eV 
is not expected to deviate from the speed of light by a factor on the order of 1CU 5 . It does not really matter in this 
case that the OPERA experiment has measured a deviation of v v from c in the superluminal direction. A hypothetical 
experimental result for v v — c < in the sw&luminal direction, of the same order-of-magnitude, as indicated in Eq. (14), 
would have been equally surprising. According to previous neutrino data [qTU"2"]. OPERA was not expected to find a 
deviation \v v — c| in the neutrino propagation velocity of the order-of-magnitude given in Eq. ( 15 ) . In light of Eqs. \12\ 



and (13), the OPERA signal would otherwise correspond to a particle with a rest mass in the range of a hundred 
MeV, or, with an effective mass of the neutrino that grows linearly with the energy. Unfortunately, neither the Higgs 
mechanism nor the Gross-Neveu model, induce a mass that depends on the energy. Once the vacuum expectation 
value of the background field that generates the mass is fixed, the mass of the constituent particle is also fixed. We 
find that it is indicated to investigate genuine neutrino mass running due to interactions which have hitherto not been 
introduced into the standard model. In Sec. IV of this paper, we write down a chiral Yukawa interaction which might 
induce a neutrino mass running with the experimentally observed parameters. 



Acknowledgments 

Helpful conversations with B. J. Wundt are gratefully acknowledged. The author acknowledges support from the 
National Science Foundation and by a Precision Measurement Grant from the National Institute of Standards and 
Technology. 



E. Recami, J. Phys. Conf. Ser. 196, 012020 (2009). 

V. L. Dadykin, G. T. Zatsepin, V. B. Karchagin, P. V. Korchagin, S. A. Mal'gin, O. G. Ryazhskaya, V. G. Ryasnyi, 
V. P. Talochkin, F. F. Khalchukov, V. F. Yakushev, M. Aglietta, G. Badino, G. Bologna, C. Castagnoli, A. Castellina, W. 
Fulgione, P. Galeotti, O. Saavedra, J. Trincero, and S. Vernetto, JETP 45, 593 (1987). 
P. Adamson et al, Phys. Rev. D 76, 072005 (2007). 

G. R. Kalbfleisch, N. Baggett, E. C. Fowler, and J. Alspector, Phys. Rev. Lett. 43, 1361 (1979). 

M. Antonello et al., ICARUS Collaboration, Measurement of the neutrino velocity with the ICARUS detector at the CNGS 
beam, e-print arXiv:1203.3433vl. 

R. G. H. Robertson, T. J. Bowles, G. J. Stephenson, D. L. Wark, J. F. Wilkerson, and D. A. Knapp, Phys. Rev. Lett. 67, 
957 (1991). 

K. Assamagan, C. Bronnimann, M. Daum, H. Forrer, R. Frosch, P. Gheno, R. Horisberger, M. Janousch, P.-R. Kettle, T. 

Spirig, and C. Wigger, Phys. Lett. B 335, 231 (1994). 

W. Stoeffl and D. J. Decman, Phys. Rev. Lett. 75, 3237 (1995). 

K. Assamagan, C. Bronnimann, M. Daum, H. Forrer, R. Frosch, P. Gheno, R. Horisberger, M. Janousch, P. R. Kettle, T. 
Spirig, and C. Wigger, Phys. Rev. D 53, 6065 (1996). 

C. Weinheimer, B. Degen, A. Bleile, J. Bonn, L. Bornschein, O. Kazachenko, A. Kovalik, and E. Otten, Phys. Lett. B 460, 
219 (1999). 

V. M. Lobashev, V. N. Aseev, A. I. Belesev, A. I. Berlev, E. V. Geraskin, A. A. Golubev, O. V. Kazachenko, Y. E. 
Kuznetsov, R. P. Ostroumov, L. A. Ryvkis, B. E. Stern, N. A. Titov, S. V. Zadorozhny, and Y. I. Zakharov, Phys. Lett. 
B 460, 227 (1999). 

A. I. Belesev, E. V. Geraskin, B. L. Zhuikov, S. V. Zadorozhny, O. V. Kazachenko, V. M. Kohanuk, N. A. Lihovid, V. M. 
Lobasheva, A. A. Nozik, V. I. Parfenov, A. K. Skasyrskaya, E. A. Sudachkov, N. A. Titiov, and V. G. Usanov, Phys. At. 
Nucl. 71, 449 (2008). 



see the URL http://cupp.oulu.fi/neutrino/nd-mass.html 



A. Chodos, A. I. Hauser, and V. A. Kostelecky, Phys. Lett. B 150, 431 (1985). 

T. Chang, A new Dirac-type equation for tachyonic neutrinos, e-print arXiv:hep-th/0011087 

T. Chang, Nucl. Sci. Technol. 13, 129 (2002). 

U. D. Jentschura and B. J. Wundt, Pseudo-Hermitian Quantum Dynamics of Tachyonic Spin-1/2 Particles, submitted, a 
preliminary version can be found as e-print [arXiv: 11 10.4171] 

A. Sommerfeld, Zur Elektronentheorie. HI. Ueber Lichtgeschwindigkeits- und Ueberlichtgeschwindigkeitselektronen, 
Nachr. k. Ges. Wiss. Gottingen, Heft 3, pp. 201-235 (1905). 
G. Feinberg, Phys. Rev. 159, 1089 (1967). 

O. M. P. Bilaniuk, V. K. Deshpande, and E. C. G. Sudarshan, Am. J. Phys. 30, 718 (1962). 



13 



J. Bandukwala and D. Shay, Phys. Rev. D 9, 889 (1974). 

H. van Dam, Y. J. Ng, and L. C. Biedenharn, Phys. Lett. B 158, 227 (1985). 

L. Xinzhou and L. Jizong, J. Phys. A 20, 6113 (1987). 

O. M. Bilaniuk, J. Phys. Conf. Ser. 196, 012021 (2009). 

S. K. Bose, J. Phys. Conf. Ser. 196, 012022 (2009). 

P. Defraigne and G. Petit, Metrologia 40, 184 (2003). 

J. Levine, Metrologia 45, 184 (2008). 

A. G. Cohen and S. L. Glashow, Phys. Rev. Lett. 107, 181803 (2011). 

G. Amelino-Camelia, Relativity in space-times with short- distance structure governed by an observer-independent (Planck- 
ian) length scale, e-print arXiv:gr-qc/0012051 
G. Amelino-Camelia, Symmetry 2, 230 (2010). 

G. Amelino-Camelia, G. Gubitosi, N. Loret, F. Mercati, G. Rosati, and P. Lipari, OPERA-reassessing data on the energy 
dependence of the speed of neutrinos, arXiv:1109.5172| 

G. Amelino-Camelia, L. Freidel, J. Kowalski-Glikman, and L. Smolin, OPERA neutrinos and deformed special relativity, 
larXiv:1110T052Tl 

X.-J. Bi, P.-F. Yin, Z.-H. Yu, and Q. Yuan, Phys. Rev. Lett. 107, 241802 (2011). 
R. Cowsik, S. Nussinov, and U. Sarkar, Phys. Rev. Lett. 107, 251801 (2011). 



L. Gonzalez-Mestres, Astrophysical consequences of the OPERA superlummal neutrino, arXiv:1109.6630 
M. E. Arons and E. C. G. Saudarshan, Phys. Rev. 173, 1622 (1968). 
J. Dhar and E. C. G. Saudarshan, Phys. Rev. 174, 1808 (1968). 

R. I. Sutherland and J. R. Shepanski, Phys. Rev. D 33, 2896 (1986). 



S. Mohanty and S. Rao, Constraint on super- luminal neutrinos from vacuum Cerenkov processes, |arXiv:1111.2725| 

M. Li, D. Liu, J. Meng, T. Wang, and L. Zhou, Replaying neutrino bremsstrahlung with general dispersion relations, 

larXiv:lllL329l 

U. D. Jentschura, Phys. Rev. A 64, 013403 (2001). 
W. Pauli, Rev. Mod. Phys. 15, 175 (1943). 

C. M. Bender and S. Boettcher, Phys. Rev. Lett. 80, 5243 (1998). 

C. M. Bender and G. V. Dunne, J. Math. Phys. 40, 4616 (1999). 

C. M. Bender, S. Boettcher, and P. N. Meisinger, J. Math. Phys. 40, 2201 (1999). 

C. M. Bender, D. C. Brody, and H. F. Jones, Phys. Rev. Lett. 89, 270401 (2002). 
A. Mostafazadeh, J. Math. Phys. 43, 205 (2002). 

A. Mostafazadeh, J. Math. Phys. 43, 2814 (2002). 
A. Mostafazadeh, J. Math. Phys. 43, 3944 (2002). 
A. Mostafazadeh, J. Math. Phys. 44, 974 (2003). 

D. G. Boulware, Phys. Rev. D 1, 2426 (1970). 

N. Moiseyev, Phys. Rep. 302, 211 (1998); C. W. McCurdy and C. K. Stroud, Comput. Phys. Commun. 63, 323 (1991). 
U. D. Jentschura, A. Surzhykov, M. Lubasch, and J. Zinn- Justin, J. Phys. A 41, 095302 (2008). 
C. Itzykson and J. B. Zuber, Quantum Field Theory (McGraw-Hill, New York, 1980). 
U. D. Jentschura and B. J. Wundt, Eur. Phys. J. C 72, 1894 (2012). 

V. K. Oikonomou, The 2d Gross-Neveu Model for Pseudovector Fermions and Tachyonic Mass Generation, e-print 
larXiv:1109.6"T70l 

I. M. Suslov, Zh. Eksp. Teor. Fiz. 120, 5 (2001), [JETP 93, 1 (2001)]. 

I. M. Suslov, Pis'ma v ZhETF 74, 211 (2001), [JETP Lett. 74, 191 (2001)]. 

O.-M. Bilaniuk and E. C. G. Sudarshan, Nature (London) 223, 386 (1969). 



